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MODELING SUBSTRATE NOISE COUPLING 
USING SCALABLE PARAMETERS 

5 

Cross Reference to Related Application 

This application claims the benefit of U.S. Provisional Application No. 
60/417,518, filed October 9, 2002, and U.S. Provisional Application No. 60/422,145, 
10 filed October 28, 2002, both of which are incorporated herein by reference. 

Background 

In modern integrated-circuit design, it is possible to integrate digital, analog, and 
RF circuitry onto a single chip. As a result of such integration, however, digital 

1 5 circuitry can inject noise through the common substrate that affects the sensitive analog 
and RF circuitry. This substrate noise coupling can severely degrade the performance 
of noise-sensitive circuits. Accordingly, substrate-noise-coupling analysis is becoming 
an important part of the design flow of integrated circuits. Typically, the effects of 
substrate coupling at low frequencies (e.g., up to 2-3 GHz) can be verified by 

20 computing the substrate resistances between all circuit parts that inject noise into the 
substrate or that are sensitive to noise. 

Some of the most commonly used methods to analyze substrate noise coupling 
involve costly trial-and-error procedures due to the lack of an efficient substrate- 
network extractor for practical circuits. Consequently, integrated-circuit designs are 

25 often delayed and engineering time increased. To solve this problem, other techniques 
have been proposed for computing substrate noise coupling. 

For instance, some known techniques involve a detailed numerical analysis of 
substrate noise coupling. In one approach, a device simulator is used to obtain a full 
numerical simulation of currents and potentials in the substrate. In this method, the 
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entire substrate is discretized, or meshed, into smaller data points, thereby creating large 
resistance and capacitance matrices. The mesh definition plays an important role in this 
technique because it involves a tradeoff between accuracy and computational efficiency. 
Consequently, this technique is impractical for large circuit designs. 
5 Another known method for calculating the substrate network is the so-called 

boundary-element method ("BEM"). The BEM requires meshing only for the contacts 
and results in a small, but dense, Z matrix. In order to calculate substrate resistances, 
the inverse Z matrix has to be computed. Dense matrix inversion, however, has a 
computational complexity of 0(N 3 ) where N is the matrix dimension. Consequently, 

10 this method is computationally intensive for large circuit designs. 

Another technique (the "preprocessed BEM"), involves first obtaining Z 
parameters from polynomial curve fitting. This method provides a simpler estimation of 
Z parameters, but requires models for different geometries and spacings to be computed 
and stored in a design-tool library. These libraries have to be adapted to each process. 

15 Accordingly, the preprocessed BEM is not efficient for large designs. 

An alternate approach to these methods uses a scalable resistance-based model 
for substrate-resistance extraction. In one such approach, a scalable resistance-based 
model is used to predict substrate noise coupling between contacts in a heavily doped or 
lightly doped CMOS process. This scalable model, however, was developed for two 

20 contacts using a two-port resistive model, and cannot be used for multiple contacts or 
account for three-dimensional variations in contact architecture. 



Summary 

In view of the issues and concerns described above, various substrate modeling 
25 methods and apparatus are disclosed. One of the disclosed methods comprises 
determining scalable Z parameters associated with at least two substrate contacts, 
constructing a matrix of the scalable Z parameters for the at least two substrate contacts, 
and calculating an inverse of the matrix to determine resistance values associated with 
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the at least two substrate contacts. The substrate can be a heavily doped substrate or a 
lightly doped substrate. In some embodiments, the matrix is an N x N matrix, where N 
is the number of contacts. In other embodiments, the act of determining comprises 
dividing a first contact into smaller portions (e.g., rectangular or square portions), and 
5 determining respective Z parameters between the smaller portions and a second contact. 
In still other embodiments, three scalable Z parameters are determined for a first contact 
and a second contact of the at least two contacts. For example, a first of the scalable Z 
parameters can be a ratio of an open-circuit voltage at the first contact to an input 
current at the first contact, a second of the scalable Z parameters can be a ratio of an 

10 open-circuit voltage at the second contact to an input current at the second contact, and 
a third of the scalable Z parameters can be a ratio of an open-circuit voltage at the first 
contact to a source current at the second contact. In other embodiments, at least one of 
the scalable Z parameters is a function of contact area and contact perimeter, or at least 
one of the Z parameters is a function of contact geometry and contact separation. In yet 

1 5 another embodiment, the scalable Z parameters comprise a first set of scalable Z 

parameters associated with resistances between the at least two substrate contacts and a 
groundplane, and a second set of scalable Z parameters associated with cross-coupling 
resistances between the at least two substrate contacts. The scalable Z parameters of the 
first set can be based on a first model equation and the scalable Z parameters of the 

20 second set can be based on a second model equation. For example, in one embodiment, 
the first model equation can be 

z- ! , 

K x Area -f K 2 Perimeter + K 3 
wherein Z is a ratio of an open-circuit voltage to input current for a selected contact 
with other contacts being open circuits, Area is an area of the selected contact, 
25 Perimeter is a perimeter of the selected contact, and K\ 9 K 2 , and K$ are parameters that 
are dependent on substrate properties and, in one certain embodiment, can be 
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determined by curve fitting based on a simulation. In another embodiment, in which the 
substrate is a lightly doped substrate, the first model equation is 

z- ! , 

K x Perimeter + K 2 

wherein Z is a ratio of an open-circuit voltage to an input current for a selected contact 
5 with other contacts being open circuits, Perimeter is a perimeter of the selected contact, 
and K\ and Ki are parameters that are dependent on substrate properties and, in one 
certain embodiment, can be determined by curve fitting based on a simulation. In still 
another embodiment, the second model equation for a selected pair of contacts having a 
fixed relative position y is 

10 Z = ae"*, 

wherein Z is a ratio of an open-circuit voltage at a first contact to a source current at a 
second contact, x is a separation between the first contact and the second contact, a is a 
value of Z when x is zero, and p is a parameter that is dependent on substrate properties 
and, in one certain embodiment, can be determined by curve fitting based on a 

1 5 simulation or a measurement. In one alternative of this embodiment, the first contact 
and the second contact have a same size. In yet another embodiment, the second model 
equation for a selected pair of contacts having a fixed separation x is 

Z-ay 1 + 6y + c, 

wherein Z is a ratio of an open-circuit voltage at a first contact to a source current at a 
20 second contact, y is a relative position between the first contact and the second contact, 
and a, b, and c are scalable parameters that substantially depend on contact dimensions 
and, in one certain embodiment, can be at least partially determined by curve fitting 
based on a simulation or a measurement. In one alternative of this embodiment, a size 
of the first contact is different than a size of the second contact. In still another 
25 embodiment, the second model equation for a selected pair of contacts is 

Z = [ay 2 +by + c\e- p(x - x *\ 
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wherein Z is a ratio of an open-circuit voltage at a first contact to a source current at a 
second contact, y is a relative position between the first contact and the second contact, 
a, b y and c are scalable parameters for the substrate that depend on contact dimensions, 
x is a separation between the first contact and the second contact, x a is a value of x used 
5 in determining a, 6, and c, and P is a parameter that is dependent on substrate 
properties- In a certain embodiment, at least one of the parameters a, b 9 c or p is 
determined by curve fitting based on a simulation. In one alternative of this 
embodiment, a size of the first contact is different than a size of the second contact. In 
yet another embodiment in which the substrate is a lightly doped substrate, the second 
10 model equation for a selected pair of contacts having a fixed relative position y is 

Z = aK 0 (fix), 

wherein Z is a ratio of an open-circuit voltage at a first contact to a source current at a 
second contact, Ko is a Oth-order Bessel function of the second kind, x is a separation 
between the first contact and the second contact, and a and P are parameters that are 

15 dependent on substrate properties. In another embodiment in which the substrate is 

lightly doped, the second model equation for a selected pair of contacts predicts a value 
Z as a function of a separation x between the first contact and the second contact, 
wherein Z is a ratio of an open-circuit voltage at a first contact to a source current at a 
second contact, and log(Z) has a linear behavior when x is greater than a certain value 

20 and an asymptotic-like behavior when jc is less than the certain value. 

Another disclosed method of substrate modeling comprises determining scalable 
parameters associated with at least two substrate contacts (where at least one of the 
scalable parameters is scalable with a contact perimeter), constructing a matrix of the 
scalable parameters for the at least two substrate contacts, and calculating an inverse of 

25 the matrix to determine resistance values associated with the at least two substrate 
contacts. In one embodiment, the scalable parameters are Z parameters. In another 
embodiment, at least one of the scalable parameters is scalable with a contact 



PMB/pmb 245-66956 10/09/03 



EXPRESS MAIL LABEL NO. EV339201219US 
DATE OF DEPOSIT: October 9, 2003 



separation. The scalable parameters can comprise a first set of scalable parameters 
associated with resistances between the at least two substrate contacts and a 
groundplane, and a second set of scalable parameters associated with cross-coupling 
resistances between the at least two substrate contacts. 
5 Yet another disclosed method of substrate modeling comprises determining 

scalable parameters associated with at least three substrate contacts, constructing a 
matrix of the scalable parameters representative of the at least three substrate contacts, 
and calculating resistance values associated with the at least three substrate contacts 
from the matrix. In one embodiment of this method, the scalable parameters are Z 
10 parameters. 

Methods for determining a scalable Z parameter for a contact in a substrate, 
wherein the scalable Z parameter is associated with a resistance between the contact and 
a groundplane, are also disclosed. The methods comprise modeling the Z parameter as 
a function of contact area and contact perimeter (the function comprising at least one 

1 5 coefficient that is dependent on properties of the substrate), obtaining a plurality of 
sample data points for the Z parameter in the substrate (the sample data points being 
obtained for a range of contact sizes), and determining values of the multiple 
coefficients such that the function produces a curve that fits the sample data points. In 
one embodiment, the range of contact sizes is from about 2.4 /xm to about 100 fim. In 

20 another embodiment, the contacts are square or rectangular. In still another 

embodiment, at least a portion of the sample data points are obtained from a simulation. 
In one particular embodiment, the function is 

z= ! , 

K x Area + K 2 Perimeter + K 3 

wherein Z is a ratio of an open-circuit voltage to input current for the contact with all 
25 other contacts in the substrate being open circuits, Area is the contact area, Perimeter is 
the contact perimeter, and K u K 2 , and K3 are coefficients that are dependent on the 
properties of the substrate. 
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Methods for determining a scalable Z parameter for a pair of contacts in a 
substrate, wherein the scalable Z parameter is associated with a cross-coupling 
resistance between a first contact and a second contact of the pair of contacts, are also 
disclosed. The methods comprise modeling the Z parameter as a function of a 
5 separation x between the first contact and the second contact (the function including 
multiple coefficients, wherein at least one of the multiple coefficients is dependent on 
properties of the substrate), obtaining a plurality of sample data points for the Z 
parameter (the sample data points being obtained for a range of separations x between 
the first contact and the second contact), and determining values of the multiple 

10 coefficients such that the function produces a curve that fits the sample data points. In 
one embodiment, the first contact and the second contact have a same size. In another 
embodiment, the range of separations x comprises values of x substantially equal to or 
greater then 10 /xm. In yet another embodiment, the range of separations x is from 
about 10 /xm to about 120 /xm. In still another embodiment, at least a portion of the 

15 sample data points are obtained from a simulation. In one particular embodiment, the 
function is 

Z = ae- px , 

wherein Z is a ratio of an open-circuit voltage at the first contact to a source current at 
the second contact, a is a value of Z for jco, and P is a parameter that is dependent on the 
20 properties of the substrate. In an alternative of this embodiment, a can be determined 
from 

1 



a = - 



K l Area + K 2 Perimeter + K 3 9 
wherein Area is a combined contact area, Perimeter is a perimeter of the combined 
contact, and K\ 9 K 2 , and K3 are coefficients that are dependent on the properties of the 
25 substrate. In another alternative of this embodiment, K\ 9 K2, and K3 are determined by 
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curve fitting a to a plurality of data points obtained for a range of different Area and 
Perimeter values. 

Other methods for determining a scalable Z parameter for a pair of contacts in a 
substrate, wherein the scalable Z parameter is associated with a cross-coupling 
5 resistance between a first contact and a second contact of the pair of contacts, are also 
disclosed. The methods comprise modeling the Z parameter as a function of a relative 
position y between the first contact and the second contact where the first contact has a 
greater dimension than the second contact along a y axis (the function including 
multiple coefficients, wherein at least one of the multiple coefficients is scalable with 

10 dimensions of the first contact), obtaining a plurality of sample data points for the Z 
parameter (the sample data points being calculated for a range of positions y of the 
second contact relative to the first contact), and determining values of the multiple 
coefficients such that the function produces a curve that fits the sample data points. In 
one embodiment, the range of positions^ is from substantially zero to a length of the 

1 5 first contact along its y axis. In another embodiment, the plurality of data points are 
obtained for a contact having an area between about 2.4 /zm and 100 jtim. In still 
another embodiment, at least a portion of the sample data points are obtained from a 
simulation. In one particular embodiment, the function is 

Z-ay 1 + by + c, 

20 wherein Z is a ratio of an open-circuit voltage at the first contact to a source current at 
the second contact, y is a relative position between the first contact and the second 
contact, and a, b, and c are scalable coefficients for the substrate that depend on contact 
dimensions. In an alternative of this embodiment, c is found by 

c = ae" Px ° 

25 wherein a is a value of Z for jc 0 , p is a coefficient that is dependent on substrate 

properties, and x a is a separation between the first contact and the second contact. In 
another alternative of this embodiment, the pair of contacts is an original pair of 
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contacts, and a, b, and c are scaleable for a new pair of contacts by a ratio of Onew/ot, 
where a new is a value of a for the new pair of contacts and a is a value of a for the 
original pair of contacts. 

Still other methods for determining a scalable Z parameter for a pair of contacts 
5 in a substrate, wherein the Z parameter is associated with a cross-coupling resistance 
between a first contact and a second contact of the pair of contacts, are also disclosed. 
The methods comprise modeling the scalable Z parameter as a function of a separation x 
between the first contact and the second contact and as a function of a relative position;; 
between the first contact and the second contact where the first contact has a greater 

1 0 dimension than the second contact along a y axis (the function including multiple 
coefficients, wherein at least one of the multiple coefficients is scalable with 
dimensions of the first contact and at least one of the multiple coefficients is dependent 
on substrate properties), obtaining a first set of sample data points for the Z parameter 
(the first set of sample data points being obtained for a range of relative positions y of 

1 5 the second contact relative to the first contact for a fixed separation jc), obtaining a 
second set of sample data points for the Z parameter (the second set of sample data 
points being obtained for a range of separations x for a fixed relative position;; of the 
second contact), and determining values of the multiple coefficients such that the 
function produces a curve that fits the sample data points. In one particular 

20 embodiment, the function is 

Z = [ay 2 +by + c]e- fi (*- x '\ 
wherein Z is a ratio of the open-circuit voltage at the first contact to the source current 
at the second contact, y is a relative position between the first contact and the second 
contact, a, b, and c are scalable parameters for the substrate that depend on contact 
25 dimensions, x is a separation between the first contact and the second contact, jc a is a 
value of x used in determining a, b 9 and c, and P is a parameter that is dependent on the 
properties of the substrate. 
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Computer-readable medium storing computer-executable instructions for 
causing a computer system to perform any of the disclosed methods are also disclosed. 

These and other features are set forth below with reference to the accompanying 
drawings. 

Brief Description of the Drawings 

FIG. 1 is a block diagram schematically showing an exemplary substrate 
network having more than two contacts. 

FIG. 2A is a block diagram showing a resistive network for two exemplary 
contacts. 

FIG. 2B is a block diagram showing the resistive network of FIG. 2 A with three 
contacts. 

FIG. 3 is a block diagram showing a three-contact resistive network using a Z- 
parameter formulation. 

FIG. 4 is a block diagram showing the three-contact resistive network of FIG. 3 
using a 7-parameter formulation. 

FIG. 5 is a flowchart of an exemplary method of determining noise coupling in a 
substrate. 

FIG. 6 is a block diagram schematically showing a first exemplary heavily 
doped substrate. 

FIG. 7 is a graph showing simulated values of Zn as a function of a contact 
width in an exemplary heavily doped substrate. 

FIG. 8 is a graph showing simulated values of 1/Zn as a function of contact 
width for the exemplary heavily doped substrate of FIG. 7. 

FIG. 9A is a first image showing current-flow lines from a simulator for an 
exemplary heavily doped substrate having a single contact. 

FIG. 9B is a second image showing current- flow lines from a simulator for the 
heavily doped substrate of FIG. 9 A having two contacts with a separation of 5 jam. 
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FIG. 9C is a third image showing current-flow lines from a simulator for the 
heavily doped substrate of FIG. 9A having two contacts with a separation of 10 jim. 

FIG. 9D is a fourth image showing current-flow lines from a simulator for the 
heavily doped substrate of FIG. 9 A having two contacts with a separation of 40 |im. 
5 FIG. 10 is a logarithmic graph of Z\ 2 as a function of contact separation showing 

values calculated by a two-dimensional simulator and values predicted by an exemplary 
Z-parameter model. 

FIG. 1 1 is a logarithmic graph of Z\ 2 as a function of contact separation showing 
values calculated by a two-dimensional simulator and values predicted by an exemplary 
10 Z-parameter model for separations of 10 [xm and 0.5 |im. 

FIG. 12 is a block diagram schematically showing a second exemplary heavily 
doped substrate. 

FIG. 13 is a graph of 1/Zn as a function of contact width showing values 
calculated by a three-dimensional simulator and values predicted by an exemplary Z- 
1 5 parameter model. 

FIG. 14 is a bar graph comparing values of Z n predicted by a three-dimensional 
simulator and an exemplary Z-parameter model. 

FIG. 15 is a flowchart of an exemplary method for determining a scalable Z n 
parameter for a pair of contacts in a substrate. 
20 FIG. 16 is a logarithmic graph of Z n as a function of contact separation showing 

values calculated by a three-dimensional simulator and values predicted by an 
exemplary Z-parameter model for separations of 20 |im and 5 |im. 

FIG. 17 is a block diagram of a pair of contacts illustrating the relationship of a 
first contact to a second contact along x andj> axes. 
25 FIG. 18 is a graph of Zj 2 as a function of relative position^ for the contacts 

shown in FIG. 17 showing values calculated by a three-dimensional simulator and 
values predicted by an exemplary Z-parameter model. 
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FIG. 19A is a block diagram illustrating two exemplary contacts having a 
separation x of zero and a relative position of 100 |am. 

FIG. 19B is a block diagram illustrating the exemplary contacts of FIG. 19A 
with a relative position^ of 80 jim. 
5 FIG. 19C is a block diagram illustrating the exemplary contacts of FIG. 19A 

with a relative position^ of 50 jam. 

FIG. 20A is block diagram illustrating two exemplary contacts of equal size 
having a separation x of zero and a relative position^ of 100 fim. 

FIG. 20B is a block diagram illustrating the exemplary contacts of FIG. 20A 
1 0 with a relative position of 50 |im. 

FIG. 21 is a flowchart of an exemplary method for determining a scalable Z\ 2 
parameter for a pair of contacts in a substrate as a function of a separation x between the 
contacts. 

FIG. 22 is a flowchart of an exemplary method for determining a scalable Z n 
15 parameter for a pair of contacts in a substrate as a function of a relative position j; of the 
contacts. 

FIG. 23 is a logarithmic graph of Z\ 2 as a function of the separation x and the 
relative position^ showing values calculated by a three-dimensional simulator and 
values predicted by an exemplary Z-parameter model. 
20 FIG. 24 is a block diagram of an exemplary pair of contacts illustrating the 

relationship of a first contact to a second contact along x and y axes and showing that 
the first contact has a complex shape and can be divided into smaller portions. 

FIG. 25 is a logarithmic graph of Z\ 2 as a function of the separation;; for the pair 
of contacts shown in FIG. 24 at a separation x of 30 jam showing values calculated by a 
25 three-dimensional simulator and values predicted by an exemplary Z-parameter model. 

FIG. 26 is a logarithmic graph of Z\ 2 as a function of the separation^ for the pair 
of contacts shown in FIG. 24 at separations x of 10 jam, 30 |am, and 50 jim showing 
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values calculated by a three-dimensional simulator and values predicted by an 
exemplary Z-parameter model. 

FIG. 27 is a cross section of a three-dimensional plot of Z J2 values for the 
contacts shown in FIGS. 24 and 29 at various x and y positions. 
5 FIG. 28 is a two-dimensional plot along the x and y axes of the three- 

dimensional plot shown in FIG. 27. 

FIG. 29 is a block diagram of the exemplary pair of contacts from FIG. 24 
illustrating the relationship of the first contact to the second contact in various relative 
positions on the x wAy axes. 
10 FIG. 30 is a block diagram of an exemplary pair of contacts where one of the 

contacts is C-shaped. 

FIG. 31 is a logarithmic graph of Zi 2 as a function of a separation;; for the pair 
of contacts shown in FIG. 30 at separations x of 10 (am, 30 fxm, and 45 |im showing 
values calculated by a three-dimensional simulator and values predicted by an 
15 exemplary Z-parameter model. 

FIG. 32 is a logarithmic graph of Z\ 2 as a function of a separation x for the pair 
of contacts shown in FIG. 30 at a separation;; of 30 \xm showing values calculated by a 
three-dimensional simulator and values predicted by an exemplary Z-parameter model. 

FIG. 33 is a block diagram of an exemplary pair of contacts where one of the 
20 contacts is square. 

FIG. 34 is a graph of Zi 2 as a function of a separation;; for the pair of contacts 
shown in FIG. 33 at a separation x of 40 jim showing values calculated by a three- 
dimensional simulator and values predicted by an exemplary Z-parameter model. 

FIG. 35 is a cross section of a three-dimensional plot of Zj 2 values for the 
25 contacts shown in FIG. 33 at various x and y positions. 

FIG. 36A is a block diagram showing two exemplary contacts as they are 
subdivided into smaller portions according to a known method. 
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FIG. 36B is a block diagram showing the two exemplary contacts from FIG. 
36A as they are used in an exemplary Z-parameter model. 

FIG. 37A is a block diagram showing an exemplary pair of contacts used to 
compare a known method with an exemplary Z-parameter model. 

FIG. 37B is a block diagram showing three exemplary contacts used to compare 
a known method with an exemplary Z-parameter model. 

FIG. 38 is a block diagram showing three contacts discretized according to a 
known method. 

FIG. 39 is a bar graph showing a comparison of Z u measured from a test chip 
and predicted by an exemplary Z-parameter model. 

FIG. 40A is a block diagram showing a comparison of a Z n value for an 
exemplary L-shaped contact as measured in a test chip and as predicted by an 
exemplary Z-parameter model. 

FIG. 40B is a block diagram showing a comparison of a Zn value for an 
exemplary U-shaped contact as measured in a test chip and as predicted by an 
exemplary Z-parameter model. 

FIG. 40C is a block diagram showing a comparison of a Zn value for an 
exemplary square-shaped contact as measured in a test chip and as predicted by an 
exemplary Z-parameter model. 

FIG. 41 is a logarithmic graph of Z\ 2 as a function of a separation x showing 
values measured on a test chip and for values predicted by an exemplary Z-parameter 
model. 

FIG. 42 is a block diagram schematically showing three exemplary pairs of 
contacts having a constant separation x and different contact sizes. 

FIG. 43 is a graph of a versus contact dimension for the pairs of contacts shown 
in FIG. 42 showing values measured on a test chip and values predicted by an 
exemplary Z-parameter model. 
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FIG. 44 is a block diagram schematically showing an exemplary lightly doped 
substrate. 

FIG. 45A is a first image showing current-flow lines from a simulator for an 
exemplary lightly doped substrate having a channel-stop implant layer. 

FIG. 45B is a second image showing current-flow lines from a simulator for an 
exemplary light doped substrate having no channel-stop implant layer. 

FIG. 46 is a graph of Z u as a function of contact dimensions showing values 
simulated in a two-dimensional simulator and values predicted by an exemplary Z- 
parameter model. 

FIG. 47 is a logarithmic graph of Z i2 as a function of contact separation showing 
values simulated in a two-dimensional simulator and values predicted by an exemplary 
Z-parameter model. 

FIG. 48 is a bar graph illustrating the dependence of Z u on die dimension in an 
exemplary lightly doped substrate. 

FIG. 49 is a bar graph illustrating the dependence of Z\ \ on die dimension in an 
exemplary heavily doped substrate. 

FIG. 50 is a graph of Z\\ as a function of contact dimensions showing values 
simulated in a three-dimensional simulator and values predicted by an exemplary Z- 
parameter model. 

FIG. 51 is a logarithmic graph of Z\ 2 as a function of contact separation showing 
values simulated in a three-dimensional simulator. 

FIG. 52 is a graph of Z\ 2 as a function of contact separation showing values 
simulated in a two-dimensional simulator and values predicted by an exemplary Z- 
parameter model. 

Detailed Description 

Disclosed below are representative methods and apparatus for modeling 
substrate noise coupling. The disclosed methods should not be construed as limiting in 
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any way. Instead, the present disclosure is directed toward novel and nonobvious 
features and aspects of the various disclosed embodiments, alone and in various 
combinations and subcombinations with one another. The methods are not limited to 
any specific aspects or features, or combinations thereof, nor do the methods require 
5 that any one or more specific advantages be present or problems be solved. 

Although the operations of the disclosed methods are described in a particular, 
sequential order for convenient presentation, it should be understood that this manner of 
description encompasses rearrangement, unless a particular ordering is required by 
specific language set forth below. For example, operations described sequentially may 

10 in some cases be rearranged or performed concurrently. Moreover, for the sake of 
simplicity, the disclosed flow charts typically do not show the various ways in which 
the disclosed methods can be used in conjunction with other methods. Additionally, the 
detailed description sometimes uses terms like "determine" and "provide" to describe 
the disclosed methods. These terms are high-level abstractions of the actual operations 

15 that are performed. The actual operations that correspond to these terms will vary 
depending on the particular implementation and are readily discernible by one of 
ordinary skill in the art. 

As more fully described below, the disclosed embodiments can be applied to a 
wide variety of fabrication processes, including CMOS processes. Further, some of the 

20 methods described herein can be implemented in software, stored on a computer- 
readable medium, and executed on a computer. Some of the disclosed methods, for 
example, can be implemented in an electronic-design-automation ("EDA") tool, such as 
a design and verification tool, which includes a user interface configured to input and 
output substrate parameters and other design information. Such methods can be 

25 executed on a single computer or on a networked computer (e.g., via the Internet, a 
wide-area network, a local-area network, a client-server network, or other such 
network). For clarity, only selected aspects of software-based implementations are 
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described. Other details that are well known in the art are omitted. For the same 
reason, computer hardware is not described in further detail. 

General Description 

FIG. 1 is a block diagram showing an exemplary substrate 100 having multiple 
contacts 110. As can be seen in FIG. 1, the overall substrate-coupling-resistance 
network is a complex one. A contact 1 10 has an effective resistance 130 to a backplane 
conductor 120 as well as an effective cross-coupling resistance 132 to another contact 
111. Two-port resistance-based models are generally unable to calculate resistance 
networks such as that of FIG. 1. Such models typically do not account for resistance 
changes associated with more than two contacts. As an example, consider a two-port 
resistive-network formulation derived for a 3-contact case. FIG. 2A shows simulated 
resistance values for two identical contacts 210, 21 1 at a separation of 10 jam. As 
shown in FIG. 2B, the addition of an identical third 212 contact at a 10 jim separation 
from the contacts 210, 21 1 alters all the resistance values of the network. According to 
the two-port resistance-based models, the resistance values in FIG. 2B should be 
identical to the values in FIG. 2A. Therefore, the two-port resistance model is 
inadequate for networks having more than two contacts. 

To overcome this and other deficiencies, certain embodiments of the disclosed 
methods use scalable Z parameters that relate a voltage applied to a first contact and an 
associated current at a second contact. In general, for an arbitrary pair of contacts j and 
w, an associated Z parameter, denoted as "Zj m ," is a ratio of an open-circuit voltage at 
the contact j to the source current at the contact m. In other words, the Z parameter 
relating contacts j and m is the open-circuit voltage at contact j (Fj) divided by a source 
current at contact m (I m ) 9 or, 
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Generally, the parameter Z, m does not change due to other contacts when the separation 
between the additional contacts is greater than a certain value. For example, in certain 
substrates considered herein, the parameter Zj m does not change if additional contacts 
are separated by at least 10 jam. A contact j also has a parameter Zjj associated with a 
5 groundplane. The parameter Zjj can be defined as a ratio of open-circuit voltage at the 
contact j to a current at the contact j\ with all other contacts considered as open circuits. 

Although 7 parameters are related to Z parameters, Y parameters are unsuitable 
for describing the substrate network. Y parameters are short-circuit parameters. Thus, 
for every 7-parameter measurement, the substrate network is altered because one or 

10 more circuit nodes are grounded. Moreover, a /-parameter formulation does not take 
into account cross-coupling resistances between two grounded nodes, and resistances 
from grounded nodes to the backplane. FIGS. 3 and 4 illustrate the difference between 
the resistive networks for Z-parameter and 7-parameter formulations for a four-contact 
example. In FIG. 3, a Z-parameter formulation 310 is shown for an exemplary four- 

1 5 contact network. As can be seen, the Z-parameter formulation 310 includes resistances 
from each contact to the backplane and cross-coupling resistances between each pair. 
In FIG. 4, a 7-parameter formulation 410 is shown for the same network. In FIG. 4, 
resistors R 2 2, R33, R44, R23 and R 34 are eliminated as they are interconnected grounded 
nodes. 

20 FIG. 5 is a flowchart 500 showing one representative embodiment of a general 

method for calculating substrate noise coupling in a network having two or more 
contacts. In process block 510, scalable parameters (e.g., Z parameters) are determined 
for two or more contacts in the substrate. For instance, in one exemplary embodiment, 
Z parameters are calculated for the multiple contacts in the network. In process block 

25 520, a matrix of the scalable parameters is constructed. In one particular embodiment, 
the matrix is an N x N matrix, where N is a number of contacts in the network. For 
example, matrix entries corresponding to two contacts j and m can be included in such 
an N* N matrix as: 
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• H 


• " Z jm * • • 




• • z mm * • • 



In process block 530, resistance values are determined from the scalable 
parameter matrix. In one embodiment, for instance, the inverse of the matrix is 
calculated and used to determine the resistance values. The resulting resistance values 
5 are indicative of substrate noise coupling and can be evaluated by a designer or design 
tool to determine whether any modifications should be made to the design. 

Modeling Substrate Noise Coupling in Heavily Doped Substrates 

In order to explain how multiple, three-dimensional geometries influence the 

10 calculation of scalable parameters in a heavily doped substrates, an examination of 
representative Z-parameter models in two dimensions is helpful. Two-dimensional 
models for substrate noise coupling based on Z parameters can be developed using two- 
dimensional device simulators, such as MEDICI, Version 2000.2.0, available from 
Synopsys Corporation. In general, the use of a simulator provides accurate results, 

15 allows a user to define nonhomogenous substrate layers, and provides insight into 

substrate-coupling mechanisms for different fabrication processes. As described above, 
simulators are computationally intensive and impractical for large designs. Two- 
dimensional device simulators typically use process information, such as doping 
concentrations, layer thicknesses, and the size and separation of the injecting and 

20 sensing contacts. Further, certain two-dimensional simulators (e.g., MEDICI) can 
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generate current-flow lines, equipotential lines, and Y or Z parameters from the process 
information. 

For purposes of illustrating a representative two-dimensional Z-parameter 
model, reference is made to an exemplary heavily doped substrate process. This 
5 substrate process is used for illustrative purposes only and should not be construed as 
limiting the disclosed technology. Instead, the disclosed principles and methods can be 
applied to any substrate process. For example, as more fully explained below, the 
disclosed principles and methods can be extended to lightly doped substrates. FIG. 6 
shows a heavily doped substrate 600 comprising a heavily doped p + channel-stop 

10 implant layer 610, a lightly doped epi layer 612, and a heavily doped p + bulk region 
614. The resistivities and thicknesses for these layers are shown in FIG. 6. For 
purposes of this disclosure, the symbol "w" is sometimes used to indicate micrometers. 

FIG. 7 shows a graph 700 of Z u as a function of contact width w\ in a range of 
0.5 jam to 80 ^m. FIG. 8 shows a graph 800 of 1/Zn as a function of w\. As can be 

15 seen by curve 710 in FIG. 7 (based on data from a two-dimensional simulation), an 
exponential relationship exists between Z\\ and contact width. Similarly, curve 810 in 
FIG. 8 shows a linear relationship between 1/Zn and the contact width. According to 
one exemplary method, this relationship is modeled using the following the linear 
equation: 

20 — = a+bw [9 (3) 

wherein a and b are process parameters that can be obtained from curve fitting, and w\ 
is the contact width. For example, for the exemplary substrate profile shown in FIG. 6, 
the following values were obtained: 

a = 6.5942 x 10~ 5 — b = 3.5092 — . (4) 
Q Qm 

25 Data points 812 in FIG. 8 show simulated values of 1/Zn at selected widths w\ 9 whereas 
curve 810 shows predicted values of 1/Zn for an arbitrary width w\. As can be seen 
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from FIG. 8, the model of Equation (3) closely approximates the data extracted from the 
two-dimensional simulation. 

Another aspect of modeling Z\\ is to determine the influence of neighboring 
contacts on Zn. FIGS. 9A-9D show current-flow lines for different separations between 
5 two contacts in an exemplary substrate substantially similar to the substrate shown in 
FIG. 6. The current-flow lines shown in FIGS. 9A-9D provide insight into the 
relationship between Z u and the contact locations. In particular, FIG. 9A shows two- 
dimensional simulation results for an exemplary single contact, FIG. 9B shows two 
exemplary contacts at a 5 jam separation, FIG. 9C shows the contacts at a 10 |um 

10 separation, and FIG. 9D shows the contacts at a 40 jam separation. For each of these 
simulations, injector-contact size was w\ = 0.5 ^m and sensor-contact size as wi = 10 
^im. From the two-dimensional simulation, the following Z parameter values were 
obtained: for FIG. 9A, Z n = 14.79 KG; for FIG. 9B,Z U = 13.05 Kfi; for FIG. 9C, Z n = 
14.18 K«; and for FIG. 9D, Z n = 14.79 KQ. As can be seen in FIGS. 9A-9D, current 

15 distribution is a function of sensor-contact separation. Consequently, values of Z\\ also 
depend on separation. In FIG. 9D, the current-flow lines (current distribution) are 
substantially the same as in the single-contact case. Thus, for separations larger than a 
certain value (in the illustrated substrate, about 10 jxm), the injector contact behaves as a 
single contact, and Z\ \ is unaffected by a neighboring contact. For instance, in this 

20 example, the value of Z u for a selected contact is within about 1% of its isolated value 
for 20 jxm or greater separation from other contacts and is within 10% of the isolated 
value for separations of 10 |im or more from other contacts. Accordingly, in one 
exemplary embodiment, the Zu model described above is used for contacts separated 
from other contacts by about 10 jam or larger. In other embodiments, however, the 

25 separation for which this Z\ \ model is used can be different (because of different 
substrate characteristics, a different desired accuracy, or for other reasons). 
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10 



15 



FIG. 10 shows a logarithmic graph 1000 for values of Z\ 2 for various separations 
between source and sensor contacts. In particular, graph 1000 shows data points 1010 
obtained from a two-dimensional simulator. From FIG. 10, it can be seen that Z\ 2 is an 
exponential function of separation that can be modeled as, for example: 

Z l2 =ae-* (5) 
wherein a and P are process parameters obtained from curve fitting. Curve 1012 in 
FIG. 10 was plotted according to Equation (5) and shows good agreement with data 
points 1010. 

FIG. 1 1 is a logarithmic graph 1 100 of Z J2 as a function of contact separation. 
In particular, graph 1 100 shows data points 1 1 10, 1 1 12 that correspond to simulation 
results for contact widths of 10 jam and 0.5 jam, respectively, and curves 1 120, 1 122 
that correspond to values predicted by Equation (5). From FIG. 1 1, it can be seen that p 
(the slope of curves 1 120, 1 122) is substantially independent of the contact widths. For 
example, the values of a and p obtained for the exemplary heavily doped substrate 
shown in FIG. 6 are: 



a = 233 Q £ = 1.0666xl0 5 — . 

m 



(6) 



The cross-coupling resistance values predicted by the disclosed Z-parameter 
model can be verified by comparing the resistance values obtained from the model with 
simulated values. The two-port Y parameters for the substrate model are given by: 



20 


Y = 






i 


Z 22 


-Z, 2 




G„ +G n 










~|z| 


_- Z 21 






-G 21 



-G n j 

^21 + ^22 J 



(7) 



Therefore, 



(8) 



which gives, 
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^.2=^ 1 = ^r L (9) 

Z 12 z 21 

wherein \y\ is a determinant of the Z matrix. The 7, Z, G and matrices are symmetric 

because the resistive network between the contacts is a reciprocal network. 

Table 1 below shows the R u values for different separations (x) and contact 
5 widths calculated using the Z-parameter model described above and using a two- 
dimensional simulator. 



Source 
Width (urn) 


Sensor 
Width (urn) 


Ri 2 (x=10um) 


R12 (x=40 urn) 


R12 (x=100 jam) 


Model 


Sim. 


Model 


Sim. 


Model 


Sim. 


0.5 


0.5 


39.2KQ 


38.7KQ 


1.08MQ 


1.04MQ 


659MA 


591MQ 


0.5 


5 


40KQ 


38.6KQ 


1.08MQ 


1.04MQ 


651 MO 


590MQ 


0.5 


100 


43.1KO 


38.6KO 


1.08MQ 


1.04MQ 


649M0 


590MO 


2 


5 


40.2KQ 


38.5KQ 


1.08MO 


1.04MQ 


649MQ 


589MO 


5 


5 


40.6KQ 


38.5KG 


1.08MO 


1.04MQ 


647MQ 


588MQ 


10 


100 


41.9KQ 


38.5KQ 


1.04MQ 


1.04MQ 


628MQ 


588MQ 



Table 1: Comparison of Z-parameter model and simulation R p values 

10 

Table 1 shows that the largest difference between the resistance values predicted by the 
Z-parameter-based model and the simulated results is 1 1 .66% for a separation of 10 jam, 
while for the 40 jxm and 100 |im separations, the errors are 3.8% and 1 1 .5%, 
15 respectively. 

Although the two-dimensional model is useful to illustrate the relevant 
substrate-coupling mechanisms, it cannot account for the influence of three-dimensional 
geometries on substrate noise coupling. A scalable, three-dimensional, Z-parameter 
model can be produced using a three-dimensional simulator {e.g., MOMENTUM EM, 
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10 



Advanced Design System 1.5, from Agilent Technologies). A typical three- 
dimensional simulator allows a user to define different contact shapes and produces an 
Alport, S-parameter output file, wherein AT is a number of contacts in the layout. 

Z parameters can be calculated from S parameters in the following manner. 
The S and Y parameters for a 2-port network are given by: 

bx^Snax+Snai h = Y n V x + Y n V 2 (10) 

bi = £1201 + $2202 * 2 = ^12^1 + ^22^2 

Substituting for a { j and by in terms of the port voltages and currents in the S-parameter 
equation: 



and 



V 2 - Z 0 i 2 _ e 

2j& ~ S2 ' 



f \ 

Vi+Z 0 i x 



2JZ, 



+ S 



22 



0 J 



V 2 +Z 0 i 2 ^ 



2JZ, 



0 ) 



From Equation (11): 



(11) 



(12) 



Z 0 h + Si \Z4\ + Si 2 Zoi 2 = V\ -S\\V\ - S i2 V 2 . 



(13) 



1 5 From Equation (12): 



Z 0 i 2 + S 2 \Zoh + S 22 Zoi 2 =V 2 - 5 21 V\ - S 22 V 2 . 



(14) 



Then, 



S u +1 5 12 

] ^22 ^ 





i, 




1- 




.h. 







Ml 



M2 



Z 0 (S + /) 



5 2 i 1 - 5 22 



'2j 



L'2J 



= (i-s) 



iA 



= Y 



L V 2J 



(15) 
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So: 

Y = ±(S + I)-\I-S) (16) 
z o 

and 



Z=T\ (17) 

5 For the special case of two contacts, the Z and 7 parameters are related to each 

other as follows: 

7 -XlL 7 - ~ y i2 7 _ ~ Y 2\ 7 _ ill / 1 o\ 

^11 - |y| » ^12 - |y| > ^21 - |y| 5 ^22 - jpj l 15 J 

For purposes of illustrating the three-dimensional model, the substrate shown in 
FIG. 12 is considered. Substrate 1200 in FIG. 12 is a 0.5 jxm MOSIS HP CMOS 
10 process. The resistance values used for the examples discussed herein comprise 

average values for the MOSIS HP CMOS process and differ somewhat from a general 
heavily doped substrate. 

Data obtained from three-dimensional simulations (e.g., MOMENTUM) 
indicates that Z n is a function of both contact area and contact perimeter. Accordingly, 
1 5 in one exemplary embodiment, the following equation can be used to model Z\ \ : 

Z u = (19) 

K x Area + K 2 Perimeter + K 3 

where K u K 2 and AT 3 are empirically fit parameters, which are discussed more fully 
below. 

FIG. 13 shows a graph 1300 of \IZ\\ as a function of contact width for square 
20 contacts. As can be seen by curve 1320 based on Equation (19), 1/Zn increases 
quadratically with contact width. Further, various data points 1310 calculated by a 
three-dimensional simulator show that the curve 1320 predicted by Equation (19) is 
accurate. FIG. 14 shows a bar graph 1400 comparing results from the above Z- 
parameter model (labeled as "model" in the graph 1400), simulation results (labeled as 
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"simulation"), and results from measurements made in a test chip (labeled as 
"measurement") for a variety of different contact dimensions. FIG. 14 shows that the 
predicted results using Equation (19) are superior. Minor deviations between the 
modeled and simulated Z\\ values and the measured values are expected due to the 
5 uncertainty in the substrate-doping profiles used in the model and the simulations. 

The dependence of Z\\ on nearby contacts can also be examined for three- 
dimensional geometries. Three-dimensional simulations (in the illustrated graphs, from 
MOMENTUM) show that for contact separations larger than a certain value, Z\\ values 
do not change due to nearby contacts. In particular, simulation results show that for a 

10 square 5 jam contact with no surrounding contacts, the Z n value is 535 Q; for a 5 jam 
contact surrounded at each side by four identical contacts with a separation of 0.5 jam, 
the value of Z\ \ drops to 403 Q; and for a 5 jxm contact surrounded by four identical 
contacts with a separation of 10 jam, the value of Z u is 528 Q. From these and other 
simulations, it was determined that for the exemplary substrate shown in FIG. 12, the 

1 5 value of Z\ \ converges to a single contact value at a separation x substantially equal to 
or larger than 10 jam. 

FIG. 15 is a flowchart 1500 showing one exemplary method of determining a 
Z\\ parameter for a given contact. At process block 1510, the Z parameter is modeled as 
a function of the area and perimeter of the contact, where the function also includes 

20 multiple coefficients (or parameters), at least some of which are dependent on the 
properties of the given substrate. For example, the Z parameter can be modeled 
according to Equation (19) with the multiple coefficients being K u K 2 and AT 3 . At 
process block 1520, a plurality of different data points (e.g., ten) are calculated for Z n 
using a single contact (e.g., having a square and/or rectangular geometry). The data 

25 points can be obtained, for example, from a three-dimensional simulator or any other 
suitable method. In one exemplary embodiment, the contact sizes are varied over a 
range of values. For example, contact sizes can be varied from 2.4 [im to 100 Jim, 
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though this range is not limiting in any fashion. At process block 1530, the values of 
the multiple coefficients in the function (e.g., K\, K 2 and K 3 ) are determined by curve- 
fitting the function with the sampled data points. For example, in one exemplary 
embodiment, the values of K\, K 2 and K 3 are determined by curve-fitting Equation (19) 
5 to the Z\ \ data. A variety of known curve-fitting techniques can be used to obtain the 
relevant values. For the MOSIS HP CMOS process illustrated in FIG. 12, for instance, 
the following values for K u K 2 , and K 3 were obtained: K\ = 3.191 1 x 10 6 (1/Qm 2 ); K 2 = 
47.6177 (1/Qm); and K 3 = 7.0579 x 10 -4 (1/Q). 

In contrast to Z\\, the value of Z\ 2 is typically a function of the contact 

10 geometries and spacing. FIG. 16 is a logarithmic graph 1600 showing the value of Z !2 
as a function of separation*. In particular, graph 1600 shows simulated Z\ 2 values at 
data points 1610 for 5 Jim square contacts and at data points 1620 for 20 jun square 
contacts. As can be seen in FIG. 16, the values Z\ 2 have an exponentially decaying 
behavior with increasing separation x. According to one exemplary embodiment, this 

1 5 behavior can be modeled as: 

Z n =ae~ fix (20) 
wherein p is a process-dependent parameter, and a is a value of Z\ 2 at x = 0. The value 
of Of generally depends on contact dimensions. At zero separation, the two contacts 
merge into a single contact. The value of a is equal to the Z u value of this single, 

20 merged contact. Therefore, the value of a can be calculated using the total area and 
perimeter of the merged contact. In this model, a introduces the area and perimeter 
dependence of contacts into the Z\ 2 model so that the contact cross-coupling resistances 
depend on the area, perimeter, and spacing of contacts. In contrast, the previous models 
completely ignore perimeter dependence for cross-coupling resistances. 

25 From simulations, it has been observed that p is independent of the contact 

dimensions and is a constant for a given substrate process. For this reason, the value of 
P can be obtained by curve-fitting the value to simulated or measured data. Predicted 
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values from Equation (20) are also shown in FIG. 16. In particular, curve 1630 shows 
the predicted values for 5 /xm contacts, whereas curve 1640 shows the predicted values 
for 20 fim contacts. As can be seen in FIG. 16, there is good agreement between the 
simulated values and the predicted values. It can also be seen that p (the slope of the 
5 curves 1630 and 1640) is independent of contact geometry since the slopes of the curves 
1630, 1640 for different contact sizes are about the same. 

The ability to model contacts having various shapes and sizes relative to one 
another is significant, as most circuits do not have uniform, simple contacts. An 
example of two contacts with different sizes is shown in FIG. 17. In particular, FIG. 17 

10 shows a first contact 1710 separated from a second contact 1720 by a separation x 9 
wherein x is the distance between inner edges of the two contacts, and a relative 
position j; that is defined to be zero when the bottom edges of two contacts are aligned, 
and increases in the direction of the arrow 1730. According to the disclosed Z- 
parameter model, Z\ 2 between two contacts decreases exponentially with an increase in 

15 the separation x, as was illustrated in FIG. 16, wherein both contacts were assumed to 
have the same length in the j> direction. 

FIG. 18 is a graph 1800 showing the value of Z\ 2 versus relative position;; for a 
fixed separation of 35 /xm and a width of 100 /xm for the contact 1710 and a width 10 
/xm for the contact 1720. As seen in FIG. 18, the value of Z n reaches a maximum at j; = 

20 45 /xm, which corresponds to a position at which the center of the smaller contact 1720 
coincides with a center of the larger contact 1710 (i.e. 9 when the center coincides with 
w/2). 

According to one exemplary embodiment, the value of Zn can be modeled as a 
polynomial function ofy according to the following equation: 

25 Z n =ay 2 +by + c (21) 

wherein coefficients a and b are related to each other since Zn is symmetrical with w/2. 
Accordingly, 
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(22) 



w _-b 
2~~ 2^~* 

In Equation (21), the value of c is equal to the value of Z\ 2 at a fixed separation x = x z . 
The value of c can be found, for example, by: 

c = Z l2 \ x=x =ae- fix *. (23) 



5 For Equation (23), x a is the separation used for curve fitting the y dependence of Z ]2 . 

Hence, only one additional parameter is required to model the Zn dependence on the 

relative contact position. 

The coefficients a, b and c are scalable with contact dimensions. Therefore, 

once the parameters are extracted for a specific contact geometry, they can be scaled for 
10 different contact geometries. For instance, FIGS. 19A-19C show two contacts 1910 and 

1920 at x a = 0 for three different values of y: 100 /xm, 80 /xm, and 50 /xm, respectively. 

In FIGS. 19A-19C, the first contact 1910 has dimensions of 10 /xm x 1 10 /xm and the 

second contact 1920 has dimensions of 10 /xm x 10 /xm. For these three different cases, 

the maximum deviation in the simulated value of Z\\ is 1 .6 %. Thus, the value of Z u is 
15 substantially constant for contacts having a fixed area and perimeter. Because the 

Zn(a) value is independent ofy, coefficients a, 6, and c are scalable. Further, FIGS. 

19A-19C show that the coefficients a, b 9 and c each have the same dependence on 

contact dimensions. Thus, the coefficients can be scaled for different geometries using 

the following relationships: 

20 a new = a old ^ b new = b oU x^ c new = c old x^ (24) 

a a a 

where a 0 \d and a new are equal to the Z\\ value of the merged contacts before and after 

scaling, respectively. 

In FIGS. 19A-19C, the length of the second contact 1920 is small with respect to 

the first contact 1910. By contrast, FIGS. 20A-20B show two equal-sized contacts 

25 2010, and 2020 for two different y positions. In particular, the first contact 2010 has 
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dimensions of 10 /mi * 110 /mi, and the second contact 2020 has dimensions of 1 10 /mi 
x 10 /mi and are shown at y = 100 fim and j; = 500 /mi, respectively. As can be seen 
from the Zn values displayed in FIGS. 19A-19C, the value of Z u is substantially 
constant for a fixed area and perimeter. In FIG. 20, the illustrated values have a 
5 deviation of 3.3%. FIGS. 20A-20B further verify that the coefficients a, b and c are 
scalable with contact dimensions. 

FIGS. 21 and 22 are flowcharts showing two exemplary methods for 
determining Z n as a function of x (FIG. 21) and as a function of y (FIG. 22). As noted 
above with respect to FIG. 17, Z !2 is a function of both x and y, where x is the separation 
10 between the contacts and y is the relative position of the two contacts. As also noted 
above, Z X2 is the ratio of the open-circuit voltage at the first contact to the input current 
at the second current, with all other contacts as open circuits. Thus, Zn can be written 
as: 



z 

1 2 



(25) 

7 1= 0 



15 At process block 21 10 in FIG. 21, the Z parameter is modeled as a function of 

the separation x. The function can also include multiple coefficients (or parameter), 
where at least one of the coefficients is dependent on the properties of the given 
substrate. For example, the Z parameter can be modeled according to Equation (20): 

Z l2 =ae~ fix (20) 
20 where p is constant for a given substrate and has units of (1/m), and a depends on the 
contact dimensions and has units of Q . At process block 2120, a plurality of different 
data points (e.g., ten) are calculated for Z\ 2 using two contacts (e.g., contacts having a 
square and/or rectangular geometry). The data points can be obtained, for example, 
from a three-dimensional simulator or any other suitable method. For each of the data 
25 points, the sizes of the contacts and the separations between them can be varied. For 
example, the separation can be varied over a range of values (e.g., from 10 /xm to 120 



PMB/pmb 245-66956 10/09/03 EXPRESS MAIL LABEL NO. EV339201219US 

DATE OF DEPOSIT: October 9, 2003 

-31 - 



ixm) for multiple contact sizes (e.g., in a range from 2.4 fim to 100 Jim). At process 
block 2130, the values of the multiple coefficients in the function are determined, 
wherein at least one of the coefficients is determined by curve fitting the function to the 
sampled data points. In one particular embodiment, for example, the value of a in 
5 Equation (20) is determined by solving for the value of the two contacts at zero 
separation. In particular, the value of a is equal to the Z n value of the two contacts 
when they are merged into a single contact with zero separation. Thus, the value of a 
can be obtained from Equation (19), which can be modified as follows: 

a = (26) 

K } Area + K 2 Perim + K 3 

10 where Area and Perim are the area and the perimeter of the merged contact, 

respectively. In this particular embodiment, the value of p can be determined by curve 
fitting Equation (20) above with the data points and the value of a obtained according to 
Equation (26). A variety of known curve-fitting techniques can be used. 

FIG. 22 is a flowchart showing an exemplary method for determining Z\ 2 as a 

1 5 function of y. At process block 22 1 0, the Z parameter is modeled as a function of the 
relative position^. The function can further include multiple coefficients (or 
parameters), where at least one of the coefficients is scalable with the contact 
dimensions. For example, the Z parameter can be modeled according to Equation (21): 

Z 12 =ay 2 +6y + c (21) 

20 where a, b, and c are scalable with the contact dimensions and have units of fi/m 2 , Q/m, 
and 0, respectively. At process block 2220, a plurality of different data points (e.g., 
ten) are calculated for Z n using two contacts (e.g., contacts having a square and/or 
rectangular geometry). The data points can be obtained, for example, from a three- 
dimensional simulator or any other suitable method. In one particular embodiment, the 

25 geometries of the contacts are chosen such that one of the contacts is larger than the 
other, and the value of y is adjusted for several values from 0 /mi to the size of the 
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larger contact for a fixed separation x a . This process can be repeated for contacts having 
different geometries until the desired number of data points is obtained. Moreover, the 
geometry of the contacts chosen can also vary. For example, the contact sizes can be in 
the range from 2.4 /xm to 100 /mi, though this range is not limiting. At process block 
2230, the values of the multiple coefficients in the function are determined, wherein at 
least one of the coefficients is determined by curve fitting the function to the sampled 
data points. In one particular embodiment, the value of c in Equation (21) is 
determined. The value of c can be obtained using the following relationships: 

2 



u\ y=0 =c and Z l2 \ y=Q =ay z +by + c _=c (27) 



>>=o 



10 Thus, 



c = Z 12 |^ =«-*". (28) 



The values of a and b in Equation (21) can be determined by curve fitting to the data 
points obtained in process block 2220 and the value of c. A variety of known curve- 
fitting techniques can be used. 
15 In certain embodiments, once the coefficients have been extracted for a specific 

pair of contacts, new coefficients can be obtained for any contact size. For example, 
with respect to the model shown in Equation (21), when the contact dimensions change, 
the new value of a for the contact pair can be calculated by: 

K x x Area new +K 2 x Perim new + K 3 ^ 
20 where Area new and Perim new are the area and the perimeter of the new merged contacts 
at zero separation, respectively. Likewise, Zn for different contact sizes can be 

calculated by scaling coefficients a, b and c by . 

a 

When x mdy change at the same time, Z n can be modeled as the product of 
these two effects according to the following equation: 
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Z 12 = [ay 2 + 63; + c]e"^ _Ara) . (30) 
Any of the methods described above for determining a scalable Z parameter can be 
modified for use with Equation (30). 

FIG. 23 is a logarithmic graph 2300 showing the value of Z n versus the 
5 separation x and relative position^ as they change simultaneously. In particular, graph 
2300 shows two curves 2310, 2320 representing the values of Z\ 2 predicted by a three- 
dimensional simulator and by Equation (30), respectively. Because the two curves 
2310, 2320 are in good agreement with one another, they are essentially overlapping in 
FIG. 23. 

10 FIG. 24 shows a first contact 2410 that is partially surrounded by another 

contact 2420. As shown in FIG. 24, the number of interacting {i.e., coupling) sides to 
the first contact 2410 is more than one. Consequently, the contact 2420 can be divided 
into smaller rectangular portions 2430A and 2430B. The coupling of the contact 2410 
with portions 2430A and 2430B can be calculated separately using the models 

1 5 presented above. In one embodiment, the overall Z\ 2 between the two contacts is found 
by superposing the Zn values calculated for the smaller portions. 

FIGS. 25 and 26 are logarithmic graphs 2500 and 2600 showing a comparison 
between the values of Zn and a separation j> predicted using Equation (30) and using a 
three-dimensional simulator for the contact architecture shown in FIG. 24. For both 

20 FIGS. 25 and 26, the size of contact 2410 was 10 /xm x 10 fim and the L-shaped contact 
2420 comprised two 10 /im x 100 /im rectangles. FIG. 25 shows data points 2510 
calculated by a three-dimensional simulator and a curve 2520 predicted by Equation 
(30) for the architecture shown in FIG. 24 where x 0 = 30 /mi. Similarly, FIG. 26 shows 
data points 2610, 2612, 2614 calculated by a three-dimensional simulator for x 0 = 10, 

25 30, and 50 /im, respectively. FIG. 26 also shows curves 2620, 2622, 2624 predicted by 
Equation (30) for the same values of x 0 . As can be seen from FIGS. 25 and 26, there is 
good agreement between Equation (30) and the simulated results. 
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FIG. 27 expands these results into three-dimensions and shows a cross section 
2700 of a three-dimensional plot of Z\ 2 as a function of x and;;. FIG. 28 shows the 
cross section 2700 as a two-dimensional plot 2800 along only the x andj; axes. In FIG. 
28, the value of the Z\ 2 increases as the contour lines darken. FIG. 29 is a block 
5 diagram of the contact architecture used to obtain FIGS. 27 and 28, and illustrates 
various aspects of the figures. For example, with reference to FIG. 29, consider the 
movement of the small square contact 2910 along the contour lines shown in FIG. 28. 
When the small contact 2910 is close to portion 2930A, the major component of the 
resulting coupling will be from the portion 2930A. This situation is illustrated when the 

10 square contact 2910 is located at the point labeled "A" in FIG. 29. By contrast, when 
the small contact 2910 is at the point labeled "B," the maximum coupling is due to the 
portion 2930B. At the point labeled "C," the square contact 2910 is equidistant from 
both the portions 2930A and 2930B. Thus, the coupling from both the portions 2930A, 
2930B is equally significant. Because of the contribution from both portions 2930A, 

15 2930B, the Z n value is the same as that at point A or B, even though it is not at a 
minimum distance from either portion 2930A or 2930B. As shown in FIG. 28, a 
symmetry axis 2810 exists along the diagonal of the L-shaped contact 2920 due to the 
symmetry of the contact shape. These models were also tested on asymmetric shapes 
and a good agreement was found between simulated values and values predicted by the 

20 models. 

FIG. 30 illustrates an example of three- sided coupling between a square contact 
3010 and a U-shaped contact 3020. FIG. 31 is a logarithmic graph 3100 showing a 
comparison between the value of Z X2 and a separation;; predicted by the Equation (30) 
and calculated by a three-dimensional simulator for the contact architecture shown in 
25 FIG. 30. FIG. 3 1 shows data points 3 1 1 0, 3 1 1 2, 3 1 1 4 predicted by a three-dimensional 
simulator where xo = 10, 30 and 45 /xm, respectively. FIG. 31 also shows curves 3120, 
3122, 3124 predicted by Equation (30) for the same values of x 0 . FIG. 32 is a similar 
logarithmic graph 3200 showing a comparison between the value of Z\ 2 and a 
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separation x where;; has a value of 40 /xm. FIG. 32 shows data points 3210 calculated 
by a three-dimensional simulator and a curve 3220 predicted by Equation (30) for the 
architecture shown in FIG. 30. For both FIGS. 31 and 32, the contact 3010 was 10 /xm 
x 10 /xm and the U-shaped contact 3020 comprised three 10 /xm x 100 /xm rectangles. 
5 As can be seen by FIGS. 31 and 32, there is good agreement between Equation (30) and 
the simulated results. 

FIG. 33 illustrates an example of four-sided coupling between a square contact 
3310 and a surrounding square contact 3320. FIG. 34 is a graph 3400 showing a 
comparison between the value of Z n and a separation;; predicted by Equation (30) and 

10 by a three-dimensional simulator for the contact architecture shown in FIG. 33. In 

particular, FIG. 34 shows data points 3410 that were calculated by a three-dimensional 
simulator and a curve 3420 predicted by Equation (30) for the architecture shown in 
FIG. 33. FIG. 35 is a three-dimensional plot of Zn for varying x and y values for the 
architecture shown in FIG. 33. For both FIGS. 34 and 35, the contact 3310 was 10 /xm 

15 x 10 /xm and the square-shaped contact 3320 comprised four 10 /xm x 100 /xm 
rectangles. As can be seen by FIGS. 34 and 35, there is good agreement between 
Equation (30) and the simulated results. 

Experimental Results 

20 In this section, a comparison of the accuracy of the Z-parameter models for 

heavily doped substrate with three-dimensional simulations for a wide variety of contact 
shapes and spacings is presented. These examples demonstrate that the disclosed 
embodiments are both accurate and efficient in predicting substrate resistances. In 
previous numerical-based approaches, contacts are divided into smaller panels and a 

25 large resistive network extracted. In these panel-based approaches, charge or current 
distribution is calculated for equipotential contact surfaces. Because the current 
distribution is not uniform on the surface, each contact needs to be divided into panels. 
As a result, the resistive arrays become large and require a significant computational 
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effort. By contrast, the disclosed Z-parameter models can be used directly to generate a 
compact network representation in an efficient manner. 

FIG. 36A illustrates a non-scalable approach and shows two contacts 3610, 3620 
discretized into panels. By contrast, FIG. 36B shows the same contacts 3610, 3620 as 
5 used in certain exemplary embodiments of the disclosed models. Table 2 is a 

comparison of the resistance values calculated from an exemplary embodiment of the 
disclosed models versus the resistance values calculated using simulations. As seen in 
Table 2, the error between the two less than 5%. 





Rn(0) 


R 22 (fl) 


R. 2 (kQ) 


Model 


221 


221 


21.7 


Simulation 


232 


232 


22.2 



Table 2: Comparison of simulated and modeled resistance values 
for FIGS. 36A and 36B 

The disclosed, Z-parameter models can also be used for contacts of different 
shapes. For example, consider contacts 3710 and 3720 shown in FIG. 3 7 A. Table 3 
shows the resistance values for the contacts in FIG. 37A with a separation of 40 fim as 
calculated by an exemplary embodiment of the disclosed models and by simulations. 
As seen in Table 3, the exemplary Z-parameter model accurately computes the substrate 
resistances within 5% of the simulated values. 





Ru(0) 


R 22 (0) 


Ri 2 (kfi) 


Model 


252 


294 


10.96 


Simulation 


264 


292 


11.29 



20 Table 3: Comparison of simulated and modeled resistance values for FIG. 37A 



10 



15 



PMB/pmb 245-66956 10/09/03 EXPRESS MAIL LABEL NO. EV339201219US 

DATE OF DEPOSIT: October 9, 2003 

-37- 

The disclosed, Z-parameter models can also be used for more than two contacts. 
For example, FIG. 37B shows three contacts 3730, 3740, and 3750. Table 4 shows the 
computed resistance values for the contacts in FIG. 37B. As seen in Table 4, an 
exemplary embodiment of the disclosed models accurately computes the substrates 
5 resistances within 8% of the simulated values. 



Element 


Model 


Simulation 


Ru 


249 0 


261 0 


R22 


138 0 


141 Q 


R33 


109 0 


1180 


R12 


6.9 kO 


7.3 kO 


Rl3 


21.4 kQ 


22.4 kQ 


R23 


14.4 kO 


14.2 kQ 



Table 4: Comparison of simulated and modeled resistance values for FIG. 37B 

As can be seen, the traditional approaches divide the contacts into smaller panels 
10 before extracting the resistive network for the substrate. With certain embodiments of 
the Z-parameter models, however, the resistance values can be directly extracted. For 
example, consider the network shown in FIGS. 36A and 36B. For FIG. 36A, each pair 
of contacts was discretized into a 20 x 20 matrix. The resistance values were then 
determined by inverting this dense matrix, a process that has a complexity of 20 3 if a 
15 direct method is used. Using certain embodiments of the disclosed Z-parameter models, 
however, only the inverse of a 2 x 2 matrix needs to be calculated. Next, consider the 
three contact case shown in FIG. 38. Using the known panel-based method, a first 
contact 3810 is divided into n\ panels, a second contact 3820 into n 2 panels, and a third 
contact 3830 into n 3 panels. As a result of this contact discretization, the Z matrix will 
20 be a dense matrix of size (rt\ + n 2 + n 3 ) x (n\ + n 2 + 7*3). By contrast, certain 
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10 



embodiments of the disclosed Z-parameter model have a matrix size of 3 x 3 for the 
three contacts. 

In general, the Z matrix for the panel-based approach has a size Z size of: 

( k V k 



Z „,-._ — 



(31) 



where K is the number of contacts and w, is the number of panels for the ith contact. By 
contrast, the size of the matrix for certain embodiments of the disclosed scalable 
approach have a size of K x K, which is a substantially smaller matrix. On account of 
this smaller matrix size, the disclosed Z-parameter model is computationally efficient 
and can handle a large number of contacts. 

Table 5 shows a comparison of the approaches for the examples given in FIGS. 
36, 37, and 38. 



Example 


Model 


Simulation 


FIG. 36 


2x2 


880 x 880 


FIG. 37 


2x2 


395 x 395 


FIG. 38 


3x3 


1267 x 1267 



Table 5: Size of simulated and modeled Z matrices for Tables 2-4. 



1 5 The disclosed Z-parameter model has been validated on a test chip. The test 

chip was fabricated by a 0.35 fim CMOS TSMC process through MOSIS. The chip had 
several substrate test structures. The test structures that were used for validation of the 
Z-parameter model had several p + contacts of different shapes and sizes. All the 
contacts were connected to 60 fim x 60 /on DC probe pads. The measurement setup 

20 consisted of an HP 4156B semiconductor-parameter analyzer and a CASCADE probe 
station. 
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The test chip had a downbond that made an electrical connection from the 
package cavity to the bottom of the chip. This downbond was also connected to one of 
the chip's pads, making it possible to ground the backplane during the measurements by 
using a pin. The backplane of the die was contacted to a down-bonding metal plate 
5 through the conductive epoxy. The contact between the substrate and the epoxy 

behaved as a nonlinear element, which had to accounted for in the measurement results. 
The backplane nonlinearity could be characterized, for example, using the method 
described in A. Samavedam et al, "A Scalable Substrate Noise Coupling Model for 
Design of Mixed-Signal ICs," IEEE Journal ofSolid-State Circuits, vol. 35, pp. 895- 

10 903, June 2000. 

FIG. 39 is a bar graph 3900 showing measured Z u values 3910 versus predicted 
Z\\ values 3920 from an exemplary Z-parameter model for a variety of different-sized 
rectangular contacts. The curve-fitting parameters (in this case, the parameters from 
Equation (19)) obtained from the measured Z u values were: K\ = 3.9157 x 10 6 (1/Qm 2 ), 

15 K 2 = 55.3087 (1/Qm), and K 3 = 6.9400 x 10" 4 (1/0). By comparison, the curve-fitting 
parameters obtained from the simulated data and used in the Z-parameter model were: 
K x = 3.191 1 x 10 6 (1/Qm 2 ), K 2 = 47.6177 (1/Qm), and K 3 = 7.0579 x 10" 4 (I/O). 

The test chip included many different contact geometries in addition to the 
standard rectangular contacts. Measurement results show that the exemplary Z- 

20 parameter model also predicts Z\\ values accurately for the different shapes. For 

example, FIGS. 40A-40C show an exemplary L- shaped contact, a U-shaped contact, 
and a square-shaped contact, along with the resulting resistances determined by actual 
measurements and as predicted by the exemplary model. The area and perimeter for 
each contact is also shown in FIGS. 40A-40C. As can be seen by the values displayed, 

25 the exemplary model agrees with the measurements for all three contacts. 

In order to measure a variety of Zn values, an array of contacts with different 
separations in the test chip was used. For example, eight 2.4 jun x 2.4 fim contacts 
having four different separations were considered. FIG. 41 is a logarithmic graph 4100 
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displaying the value of Z x2 as a function of separation x between the contacts. In 
particular, the graph 4100 shows data points 4110 corresponding to the measured values 
and a curve 4120 corresponding to the predicted values from an exemplary Z-parameter 
model. The curve- fitting parameters {i.e., the parameters from Equation (20)) obtained 
5 from the measured Z\ 2 values were a = 898.15 0 and p = 8.676 x 10 4 (1/m). By 

comparison, the value of P obtained from curve-fitting to three-dimensional simulation 
results for the Z-parameter model was about 8.3 x 10 4 (1/m). The a value was 
calculated by considering the two contacts as having zero separation between them and 
using Equation (19) to determine the value of Z\ \ for the merged contact. The resulting 
10 a value was 824.37 Q , which is in good agreement with the measured value of 898.15 
G. 

The dependence of a on contact dimensions can also be verified through 
measurements taken on the test chip. FIG. 42 shows three different-sized pairs of 
contacts having a 10 fim separation. The values of a can be obtained from the ratio of 
Z 

15 -px^ where P = 8 * 676 x 1q4 ( 1/m )- FIG - 43 is a graP 11 4300 showing the value of a 

versus the dimensions of the contact when the separation x is zero. In particular, the 
graph 4300 shows data points 4310 obtained from measurements from the test chip and 
data points 4320 obtained using the exemplary Z-parameter model. FIG. 43 also shows 
that the value of a is equal to the value of Z\ \ when the two contacts are merged. 

20 

Modeling Substrate Noise Coupling in Lightly Doped Substrates 

In this section, the disclosed principles are adapted for use in modeling substrate 
noise coupling in a lightly doped substrate. For purposes of illustrating the disclosed Z- 
parameter models, reference is made to two exemplary lightly doped substrates. A first 
25 exemplary type of lightly doped substrate is illustrated in FIG. 44 and has two layers: a 
p + channel-stop implant 4410 and a uniform lightly doped layer 4412. The resistivities 
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and thicknesses for these layers are also shown in FIG. 44. A second exemplary type of 
lightly doped substrate is similar to the first, but does not have the p* channel-stop 
implant layer. These substrate processes are used for illustrative purposes only and 
should not be construed as limiting in any way. Instead, the disclosed principles and 
5 methods can be applied to any lightly doped substrate process. 

In order to explain how various three-dimensional geometries influence the 
calculation of Z parameters for lightly doped substrates, an examination of a two- 
dimensional model for Z parameters is helpful. FIGS. 45 A and 45B show two- 
dimensional simulated results of current flow from a single contact to a groundplane in 

10 the two exemplary substrates discussed above. In particular, FIG. 45 A shows a 

substrate 4500 having a thin p + channel-stop implant layer followed by a lightly doped 
layer, and FIG. 45B shows a substrate 4510 having only a homogenous lightly doped 
layer. (In FIGS. 45A and 45B, a substrate thickness of 20 /xm was used instead of 675 
/zm specified by the process information because of plotting limitations in the two- 

15 dimensional simulator.) 

As can been seen from FIG. 45 A, the current- flow lines indicate that the p + 
channel-stop implant layer creates a low-resistance path for the current flow. Thus, the 
injected current first spreads over the surface before it flows to the grounded backplane. 
In this type of substrate, the resistive path from the contact to the backplane is 

20 dependent on the overall chip surface, not just the contact dimensions. By contrast, as 
can be seen from FIG. 45B, the current does not spread on the surface for a 
homogeneous lightly doped substrate, and thus Z u changes with contact width. 

To further illustrate this point, FIG. 46 is a graph 4600 of Z n as a function of the 
contact width for the lightly doped substrate having a channel-stop implant layer and 

25 obtained using a two-dimensional simulator. In particular, data points 4610 show 
values of Z\ \ extracted from simulations. As can be seen from curve 4620, which is 
interpolated from the data points 4610, Z n stays constant with contact width. FIG. 47 is 
a logarithmic graph 4700 of Z\ 2 as a function of the separation x between an injecting 
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and a sensing contact. As can be seen in FIG. 47, data points 4710 decrease linearly 
with the separation x. This linear behavior indicates that, in the two-dimensional 
environment, Zn in lightly doped substrates decreases exponentially with the separation 
as it does in heavily doped substrates. Thus, the data points 4710 from FIG. 47 can be 
5 represented as an exponential function of the separation x. For example, a curve fit on 
the data from FIGS. 46 and 47 shows that Z\ 2 in lightly doped substrates can be 
modeled as: 

Z _ -fix 
12 -ae 

where a and p are process-dependent parameters that can be extracted using a 

10 simulator. FIG. 47 also shows a curve 4720, which was modeled using Equation (32), 
and which shows good agreement with the data points 4710. FIG. 47 also illustrates 
that the value of Zn changes at a small rate with separation. Consequently, the value of 
p is smaller for lightly doped substrates than for heavily doped substrates. 

As shown in FIG. 45 A, the p + channel-stop implant layer in a lightly doped 

15 substrate causes current spreading on the chip surface. Consequently, Zn can also be a 
function of the chip area in certain lightly doped substrates. Discontinuities in the 
channel-stop implant layer due to «-well structures can prevent the current flow on the 
chip surface. Thus, the effective chip area for a contact can be smaller than the overall 
chip area. FIG. 48 is a bar graph 4800 illustrating the simulated dependence of Z\ \ on 

20 die area in lightly doped substrates when the contact size is kept constant at 50 /im * 50 
fim. By contrast, FIG. 49 is a bar graph 4900 illustrating that there is no dependence of 
Z\\ on die area in heavily doped substrates. 

One of the shortcomings of a two-dimensional device simulator (e.g., MEDICI) 
is that only the width of the contacts is taken into account during simulations. In a 

25 typical two-dimensional simulator, the third dimension of the contact is assumed to be 
infinitely long for a given contact width. As was seen in FIG. 49, differences in die area 
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do not create a significant impact for heavily doped substrates. As was seen in FIG. 48, 
however, chip area can have a significant impact for lightly doped substrates. 

FIG. 50 is a graph 5000 illustrating the value of Z\ \ as a function of contact 
perimeter for a constant die area. As can be seen by curve 5020, which passes through 
simulated data points 5010, the value of Z\\ does have a dependence on contact size, a 
result that was not predicted by the two-dimensional simulators. 

A curve fit to the simulated results shows that Z\\ in a lightly doped substrate 
can be modeled as: 



Moreover, in certain embodiments, die-area and contact-area dependence can be 
incorporated into the Z\ \ model. 

Three-dimensional simulations for lightly doped substrates also show a different 
behavior for Zn than that predicted by the two-dimensional simulations. FIG. 51 is a 
logarithmic graph 5100 showing the value of Z\ 2 as a function of separation x for a 
lightly doped substrate. FIG. 51 indicates that Z\ 2 should be modeled such that log(Zi 2 ) 
is linear for larger separations and has asymptotic-like behavior for separations 
approaching zero. In one embodiment, for instance, a 0th order modified Bessel 
function of second kind, Ko(x), is used to model Z12. In this embodiment, the leading 
term in the asymptotic expansion of Kq(jc) for large x is: 



1 



(33) 



K x Perimeter + K : 




(34) 



The logarithm of Equation (34) is equal to: 




(35) 



And for large x, 



log(x) 



«1. 



(36) 



X 
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Hence, in Equation (34), log(x) can be neglected with respect to x and log( 



Jf )is 



constant term. Therefore, the logarithm of Ko(x) changes linearly with the separation 
for large separations. Based on these observations, in one exemplary embodiment, Z\ 2 
in lightly doped substrates can be modeled as: 
5 Z l2 (x) = aK 0 (fix) (37) 

where a and p are process-dependent parameters. These parameters can be extracted, 
for example, from a three-dimensional simulator or any other suitable method. FIG. 52 
is a graph 5200 showing Zn as a function of separation x for a lightly doped substrate. 
Data points 5210 were found using a three-dimensional simulator, whereas curve 5220 
10 was modeled by Equation (37). As can be seen from FIG. 52, the disclosed model is in 
good agreement with the simulated results. 

Having illustrated and described the principles of the illustrated embodiments, it 
will be apparent to those skilled in the art that the embodiments can be modified in 

15 arrangement and detail without departing from such principles. For example, the 

disclosed methods for substrate modeling are not limited to using Z parameters, but can 
be modified to use other network parameters instead {e.g., S parameters). In view of the 
many possible embodiments, it will be recognized that the illustrated embodiments 
include only examples and should not be taken as a limitation on the scope of the 

20 invention. Rather, the invention is defined by the following claims. We therefore claim 
as the invention all such embodiments that come within the scope of these claims. 



